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0\ ; Abstract 

^ ' The massive scalar field with A(/?^ interaction placed in (3 + 1) dimensional box is 

O . considered. The sizes of the box are V y. 13 {V = is the volume, T = is the 

^ \ temperature). The free energy is evaluated up to the 2nd order of PT. The averaging 

T— I ' on the vacuum fluctuations is separated from the averaging on the thermal fluctuations 

^ ■ explicitly. As result the free-energy is expressed through the scattering amplitudes. We 

^ \ find that in 3-loop approximation the expression for free energy coincides with the ansatz 

CN ' of Bernstein, Dashen, Ma suggested on the base of 5'-matrix formulation of statistical 

' mechanics. The obtained expressions are generalized for higher order of PT. 

o ■ 

i£} '■ 1 Introduction 

' The perturbation theory {PT) for quantum field models at finite temperature T is elaborated 

p • enough [1]. Its formal distinction from ordinary Feinman rules consists in that the loop sum- 

1^ , mations on discrete values of zero component of momentum appear (at T 7^ 0) instead of the 

• • \ loop integrations (at T = 0): 

>■ , 1 °° 



where for the boson fields 

Po = 27r/T, / = 0,±1,±2,... 

The loop integral is treated as an averaging on vacuum fluctuations. It is clear that at T 7^ 
both vacuum and thermal fluctuations of the field take place. So, the difference between the 
loop sum and loop integral is caused by thermal fluctuations alone. The standard task of the 
quantum field theory (QFT) at finite temperature is the calculation of radiative corrections to 
the free energy (or other thermodynamic quantities) of the relativistic ideal gas as functions of 
temperature, chemical potentials, renormalized masses, coupling parameters, etc. There are no 
principal difficulties, but the job is more cumbersome compared to usual Feinman diagrammar. 
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We want to consider the problem from some other standpoint. We shall try to separate the 
averaging on the thermal fluctuations from that of the vacuum fluctuations explicitly and to 
represent the contribution caused by interaction to the free energy as some integrals on Bose- 
Einstein distribution of renormalized pure QFT values. 

Note that probably the similar purpose was pursued in the papers [2,3] devoted to the 
elaboration of QFT at finite temperature in real time representation. But there was not achieved 
such simple and physically qtransparent result as we obtained in the framework of standard 
imaginary time representation. 

We shall work out the calculations for the A9? ^ QFT model (extraction of the loop integrals 
from the corresponding sums, renormalization) up to the 2nd order of PT, where the nontrivial 
3-loops diagrams appear. 



2 The 1st order of PT 

For the beginning consider the first radiative correction to the free energy. This exercise clarifies 
what result we would like to obtain. 

The partition function of the scalar field placed to the thermostat of the volume V at the 
temperature T — 1/P is 

where the action in Euclidian metrics is 

a;=(a;o,x), = + x^ □ = ^ + (V)^ 
Si[p\ = ^ j d^XLp^ix). 
The free energy density in the 1st order on A is 

/(r) = -^ = /(°)(r) + /W(r), 

where at V ^ 00, 

f'\T) = ^Sp ln[/3^K -□)] = ^ / d-V (TEln[/5^(m^ + P^)]) , 



P={Po,p), P^=Po + P^ 



/«(r) = 4<S/M> = ^f7:^/<''p(ri: ' 



Let us define the elementary diagrams to make the calculations more compact and transparent 
as follows 

O - / (^Eln[/5^(m^ , (1) 
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It is easy to compute the sum on po in (2) 

oo 

^E^^ = 7^ / ^^ + ^^?^— TT, a;=A^^^+P^. (3) 
t;;'p^ + a;2 27r i m2+p2 a;(e/?^-l)' V ^ 

— oo 

One can sec from this expression that the averaging on vacuum and thermal fluctuations cor- 
responding to one loop sum can be represented additively. By integrating eq.(3) on one 
obtains the temperature dependence of (1) i.e. 

f^'\T) = -lO=J^sJ d'pH^ - e-^n + a + bT, (4) 

where a and b are irrelevant constants (energy density and pressure of the vacuum) , their values 
depend on regularization scheme. Let us denote the vacuum and thermal loops as 

1 f d^p 



T=o (27r)'^ J w? +p^' 

2 r d^p 
(2^)3 J 2uj{eP^ - 1) ■ 

In these notations the first order correction to the free energy has the following diagrammatic 
representation 



/'"m = 5aD=^((jXD+2©2) + (7x2)). w 

Note, that (^3C3 ~ ^ O^^' divergent term (^X^ '^^) irrelevant because it 

does not depend on temperature. The first term in the r.h.s. of (6) is convergent and it 
can be represented as the integral on two-particles invariant phase volume with Bose-Einstein 
statistical factors 

8 ^JVJ 2! {2nf J 2a;i2a;2 {e^^^ - l){eP^^ -I)' ^ ' 

Keeping in mind that in the 1st order of PT the 2 — > 2 scattering amplitude is simply the 
coupling constant 

one can consider (7) as the thermal averaging of the 2-^2 scattering amplitude. Of course such 
interpretation looks artificial enough. In the 2nd order of PT the 2^2 scattering amplitude 
has nontrivial dependence on momenta, also 3 — > 3 amplitude appears. So, the next order of 
PT can bring more information about structure of representations of type (7) . Before the 

2nd order analyses the diagram ^® examined. It has nontrivial dependence on 
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temperature and is divergent. It is easy to see that this diagram can be removed by the mass 
renormahzation 

2 2 c- 2 

m — >• m — om^. 

Then the zero order diagram in (4) 

and for the cancellation of the divergent temperature dependent term of (6) is necessary 

i ®® - i^"-;® = 0. 

that is 

3 The 2nd order of PT 

In the 2nd order of PT the radiative correction to the free energy is 

/«(T) = --i,(s?M>, = -^PCCO + OD). w 

where the diagrams mean 

CXX3 = OxOxO' (10) 

^ (^/(n^'Pfc)'^''(Pl + P2 + P3 + P4)x 

^ frp3 y Sjpoi + Po2 + P03 + P04) \ 

V potio4 (^^ + Pi) (^^ + pi) (^^ + pi) (^^ + pI) y ' 

where 5-functions in (12) are mentioned as the Dirac one for integrals and the Kroneker symbol 
for sums. 

The separation of the vacuum and thermal loops from the diagram (10) is not difficult due 
to eq.(5), but the corresponding procedure for the diagram (12) requires some manipulations. 
Let us write the diagram (12) as a sum of four possible combinations of vacuum and thermal 

^""^^ ^ (^^^W) +4(^^^A^ +6(r^T) +4(j^?), (13) 




where the coefficients in the r.h.s. of (13) are chosen for the combinatorial reason: the number of 
possibilities to cut one, two or three lines connecting two vertices. The number of thermal loops 
in diagrams (13) will be equal to the number of statistical factors {e^'^ — 1)"^ in corresponding 
integrals. Letus denote the sum under the integral (12) by W {001,002, 003,004) 

por::^o4(Poi + ^i)---(po4 + ^4)' 
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and express the Kroneker 5-symbol through the integral 



fl3 

Jo 



Then calculating the simple sum 



^ipoxo cosh{xo - P /2)uj 



Po 



Po + uj^ 2a; • sinh(a;/3/2) ' 



we obtain for (14) 



W 



i fe=i \ PO 



^IPOXQ 
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Y[ f cosh(xo - (^/2)ujk 



Po + ^kJ i " k=i \ 2<^fe sinh (3u;k/2 



(15) 



Direct integration in r.h.s. of (15) gives 



W ^ 2 



+ 4- 



-I -1 



1)) 



+ 6- 



g/3 (0)3+^4) 



+ 



l[{2Me^-''- 

.k=l 

g/3(w2+W3+'^4) 

—(Vi + u;2 + COS + u;4 -^1 - 0^2 + <^3 + ^4 —(^i — u;2 — UJ3 + 104 ^ 



+ 4 



(16) 



Despite the function W is symmetric in its arguments by origin it is more convenient to represent 
the result of integration of (15) in asymmetric but compact form (16), accounting for the integral 
over momenta in (12) pick out the symmetric part of W only. Sorting (16) accordingly to the 
number of statistical factors rewrite W in the following way 



2a;i(e'3^i - 1) 2ui2uj2{e'^'^' - l){e^^-' - 1) 

4W3 



+ 



2uji2uj22uj3{ef^''^ - l)(e/^'^2 _ i)(e/3'^3 _ 1)' 



where 



/ 4 x-l 



(jji + (jj2 + + '^4 (27r)3 J V^^\ + p 



(17) 



n J^)^(^'oi + ---+Po4), (18) 



2(1120;^ 

^fc=2 

2 



(2vr) 



4 

n 

k=2 



UJl + UJ2 + UJ3 + (^4 

dpOk \ rf 

— 2 ]^\Pm H 



+ 



-UJl + tc;2 + t^s + '^4 



+ P04j 



Pi- 



rn 



W2 = 2i\{2uk 

+ 



u;i+ 002+003 + L04 
1 



+ 



+ 



+ 



n 



dpi 



Ok 



—L0i+UJ2+C^3 + ^4 CO1—LO2+UJ3+CO4J (27r)2 J V^^2 

X I 6{poi + P02 + P03 + P04) + S{-poi + P02 + P03 + P04) ) , 



X 



(20) 
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3 3 

+ 



-LOi + 102 + COS + LOi —UJi ljJ2 + U}2, — OJa 

1 3 



Thus corresponding to the notation (13) and using (17)-(21) we obtain 

@^ = ^ / (n ^'P*^) <^'(Pl + P2 + P3 + P4)W^0 



pf=p|=p|=-m2 



(2vr) 
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n 2 2 1 '^^(^^1 + P2 + P3 + Pi), 



(21) 



VV^ (27r)3 7 2a;i(e/5--l) 



(27r) 



(22) 



.2 2 
at = —m . 



The expression in the square brackets of r.h.s. (22) is none other than ordinary two-loops 
contribution to the inverse propagator, with the external momentum being on the mass shell: 



Further 



n 



(27r)6 7 \^J-^2LUk{e^^^-l), 
X Re 



2 2 ■ 

Pi — —m 



X 



(23) 





j- d^p 


r 1 


_(27r)4 J 


+ p2 


Km? + (p + pi +P2)^ 



+ 



(24) 



i. 2 2 2 
at Pi = P2 = —"m ■ 



Here in the square brackets of r.h.s. (24) one can recognize one-loop contribution (s- and 
M-channel terms) to the two-particle forward scattering amplitude on the mass shell 

(25) 

'Pi 

Just the second (-u-channel) term testifies zero angle scattering. Note that ^-channel term at 
zero angle coincides with the diagram (11) in our notations 



P2 P2 P2 

Pi pi pi 




P2~^P2 



5. -<2>- 



(26) 



Pi -^Pi 

We shall find it in the diagram (10). For the last term of r.h.s. (13) we have 



@) = (2.) 



d^Pk 



m? + (pi +P2 +P3)' 



+ 



+ 



+ {-pi+p2+P3y 



J. 2 2 2 2 

, at Pi = P2 = P3 = -m 



(27) 
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It is easy to see allowing for the momenta permutation symmetry that here the square brackets 
expression represents the contribution of four diagrams to the 3^3 zero-angle scattering 
amplitude 



P3 
P2 
Pi 



P3 

P2+P^ 
Pi 



P3 
+ Pl 
P2 



Pi 
+ P2 
P3 



P3 
P2 

Pi 



(28) 



So, we make shure that the 3-loop contribution to the free energy also exhibits the structure of 
Bose-Einstein averaging of the scattering amplitudes. 



4 Renormalization 

At the same time we see that the temperature dependent expressions (22), (24) are divergent 
owing to divergences of the vacuum loops in (23), (25), (26). Show how these divergences are 
absorbed due to mass and coupling constant renormalization 

^ - Smj - Sml X^X + SX. (29) 

The first order correction Sml is already fixed by (8). For 5ml P^^ 

SX = alA^ Sml = «2A^ (30) 

This leads to 

0-0+i©0+Y©3®+A.o. 
(X) - CX3 + ^ (rXX) + (X) . 

and for the second radiative correction (9) to the free energy we find 

- 6 ©OO - 6«i CXD + 24«20 )■ (31) 

Writing separately the contributions to f^'^\T) corresponding to different numbers of thermal 
loops we obtain: 

3 T-loops ^-yt~~^ 

^3 = 30©© +4(t^. (32) 

2 T-loops 



^2 = 9©®© -6ai©© +6((t^ - ©©©). (33) 
The expression in brackets in r.h.s. (33) is finite 

ff^ - GY^Y^ = ^'P^^'P^ X 

1 j d^p / 1 ^ 1 2 \1 

It^Y J m^ -\- p^\m^ -\- ip -\- Pi -\- P2Y m^ + (p — pi + P2Y m^ + p^J ' 



xRe 
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therefore g2 will be finite under condition of cancellation first two terms in r.h.s. (33) that is 

«i = ^<^ (35) 

With allowance for (35) we have for 

5i = 4(V^ -18 000 +24a2{T). (36) 

Using (22), (23) for the diagram (^^^ one can see from (36) that gi vanishes if 



do > - „ - , , . ^ 

4 v_yv_y 5 p ^ -m 

In usual diagram language [4] the obtained mass and coupling constant renormalization (8), 
(29), (30), (35), (37) looks as follows 

^ 2 4 6 v_/ p = — m ' 

= — ^ atpi = 0. (38) 

One can see from (38) that the renormalization obtained here from the condition of convergency 
temperature dependent part of free energy is just the same that follows from the condition of 
IPI Green's functions be finite 

^r(2)(p2) = o at p2 ^ -m^ 
r^^\pi) = -\ at pi = 0. 

Finally accounting for (31-34) the renormalized free energy in the 2nd order of PT is 
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5 Thermal averaging 



^'^3000 +4(?^ +6((^ -©©©) • (39) 



Let the contributions corresponding to averaging on 2- and 3-particles phase volume be denoted 
through /2(T) and fziT). Then from (39) we obtain 

/3(^) = (3 ©0© + 4 ) , (41) 
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and with allowance for the analytic expressions of diagrams (5), (7), (27), (34) we have for (40), 
(41) 



f2{T) = 



1 1 



2! (27r)fi 
1 1 



Re[A2-.2{pi,P2;Pi,P2)], (42) 



Re (Pi , P2 , P3 ; Pi , P2 , P3 ) ] ■ (43) 



3! (27r)9 

Some comments are to be made concerning the representation of fsit). As it follows from (27), 



(28) the diagram [T^Ty contains only 4 from 10 possible tree diagrams which saturate the 
3^3 amplitude in the 2nd order of PT. The other 6 diagrams, i.e. 



2' 
3' 

r 



3' 
2' 
1' 



3' 
1' 
2' 



V 
3' 
2' 



V 
2' 
3' 



2' 
1' 



^^44) 



have the unpleasant peculiarity to go to infinity for the forward scattering. For instance the 
first diagram from (44) is 



m'^ + (Pi +P2- P2Y {P2 - P'2Y + 2(P2 - P2)Pl 



>oo, when p2 



P'2- 



(45) 



Nevertheless the expression for (rJXEXz^) diagrams (44) give the contribution 

is finite indeed: because of 



^ dw? €) (27r)3 / 



one has 



dw? 
4 



(27r)3y 2uj{eP^-l) p^' 



d^Pk 



1 



1 



Pi. 



The way out of situation is to regularize the diagrams (44). The recipe of appropriate regular- 
ization proceeds from comparison (43), (45) and (46). Namely "zero angle" is to be understood 
as the limit 

p' = \\mi'p + e), e = {eo,e). 

The integration on phase volume is to be carried out at £ 7^ and put e = in the final result. 
Really, it is not hard to verify that in such limit the integration of the first pare of diagrams 
(44) gives 



lim 

£-♦0 



lim 

£0-^0 



f d^pi 
J 2a;i(e/5^i - 1) 
f d^pi 1 
J 2uji{eP^^ - 1) Pi 



+ 



) \w? + (pi + eY w? + (pi — e) 



in agreements with (46). 

The problem with the diagrams (44) is particular case of more general problem concerning 
n — > n amplitudes which are built up from rescattered blocks, such as 

2' 3' ••• n' 

(47) 




The diagrams (47) is singular in physical region when incoming momenta are equal to out- 
coming. The problem was intensively attacked [5,6] and regularization recipes were proposed. 
As was shown PT analyses gives the simple and natural way of solving this problem. It must 
be noted also that due to the special structure of diagrams including rescattered blocks the 
summation of subsets of such diagrams (particular case is subset of so-called ring diagrams [1]) 
removes zero angle divergences at all. 

6 Conclusion 

The remarkable simplicity and regularity of the representation (42), (43) for the contribution 
of interaction to the free energy in the 2nd order of PT gives the base for natural generalization 
to higher order of PT, namely 

oo 

/(T) - /ideal (T) = E/nm, (48) 
n=2 

where 

1 1 f ( ^ t^'^P \ 
^"^^^ ^-^.J2^J (n 2Me^^L 1) j ■■■^Pn-^P^^--- ^Pn)]- (49) 

Note that the representation for the temperature dependent part of IPI F*^^'' -function of similar 
structure was recently obtained [7]. 

In the representation (48), (49) is not hard to recognize the so-called diagrams of the first 
type in the framework of S'-matrix formulation of statistical mechanics [8]. As for the so- 
called second type diagrams we can say noting because they are constructed from bilinear 
combinations of 1mA and KeA. Corresponding terms may appear in higher (at least 3rd) 
orders of PT. 

The authors of [8] were generalized the Beth-Uhlenbeck formula and derived in the frame- 
work of nonrelativistic quantum mechanics the complete virial expansion. Being guided by 
invariant form of the representation they have supposed its validity for the relativistic case. 
Our calculations in the framework A(/?^ QFT model confirm their hypothesis with allowance 
some corrections like (n!)~^ factors in phase volume integrals reflecting the particles identity. 
Also an ambiguity in regularization procedure for the zero angle rescattering amplitudes is 
removed. 

In our opinion the (S'-matrix formulation of statistical mechanics would be good tool for the 
thermodynamic analysis of hadron matter [9-12] due to huge amount experimental information 
and well elaborated phenomenological models for the scattering amplitudes. So we consider it 
important to verify the (S'-matrix formulation of statistical mechanics by QFT methods. 

There need be no doubt that a representation analogous to (48), (49) can be derived for 
other rcnormalizable (at least normal, without confinement) QFTs because the mentioned 
representation is rather conditioned by general structure of PT diagrams than by peculiarities 
of interaction. More over, proceeding from that the thermodynamic values in this approach 
are expressed through the functionals of (S'-matrix elements one can believe it being correct 
also for hadron physics. A simple but forcible argument for this follows from the fact that 
Tr(e''^^) is invariant with respect to basis of state, so one can use in particular the set of 
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asymptotic, i.e. hadronic states. Its completeness is guaranteed by absence of any other except 
for hadrons observable asymptotic states. Unfortunately the direct derivation of formulas like 
(42), (43) from QCD, for example, is not possible because of the asymptotic states there (and 
consequently the scattering amplitudes) can not be perturbatively defined. 

The authors thanks to L.L. Jenkovszky and E.S.Martynov for fruitful discussions. 
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